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Discovery of an optical bow-shock around pulsar B0740— 28* 
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Abstract. We report the discovery of a faint Ha pulsar wind nebula (PWN) powered by the radio pulsar B0740— 28. 
The characteristic bow-shock morphology of the PWN implies a direction of motion consistent with the previously 
measured velocity vector for the pulsar. The PWN has a flux density more than an order of magnitude lower 
than for the PWNe seen around other pulsars, but, for a distance 2 kpc, it is consistent with propagation through 
a medium of atomic density nrr ~ 0.25 cm -3 , and neutral fraction of 1%. The morphology of the PWN in the 
area close to the pulsar is distinct from that in downstream regions, as is also seen for the PWN powered by 
PSR B2224+65. In particular, the PWN associated with PSR B0740-28 appears to close at its rear, suggesting 
that the pulsar has recently passed through a transition from low density to high density ambient gas. The faintness 
of this source underscores that deep searches are needed to find further examples of optical pulsar nebulae. 
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1. Introduction 



Table 1. Properties of PSR B0740-28 



Pulsar wind nebulae (PWN) around high- velocity pulsars 
provide a primary insight into the interaction between a 
pulsar and its environment. Specifically, optical observa- 
tions of such nebulae provide important information on 
pulsar velocities, and on the density, temperature and 
composition of the ambient medium. However only three 
pulsars are known to power optical bow-shock PWNe: 
B2224+65 (the 'Guitar Nebula', Cordes, et al. |1993| ), and 
the two millisecond pulsars (MSPs) B1957+20 (Kulkarni 
& Hester |l98S| ) and J0437-4715 (Bell et al. |l99§). All 
three of these pulsars have high spin-down luminosities 
and/or high space velocities. However, these pulsars dif- 
fer markedly in their spin-periods, ages and magnetic field 
strengths, highlighting the variety of pulsar winds which 
can be probed by these sources. The nebula associated 
with the neutron star RX J1856.5— 3754 further exem- 
plifies the variety of neutron stars know n to power such 
nebulae (van Kerkwijk & Kulkarni 2001 ). We have there- 
fore initiated a search for optical bow-shocks around other 
pulsars in order to characterise the properties of the asso- 
ciated pulsars, pulsar winds and ambient environments. 
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* Based on observations made with ESO Telescopes at the 
La Silla Observatory (Programme 66.D-0046) and under the 
AAT Service Mode Programme 



Spin Period" 


P (s) 


0.1667 


Period Derivative" 


P 


1.68xl0" 16 


Spin-down Luminosity" 


E (ergs s _1 ) 


1.4xl0 35 


Dispersion Measure Distance 6 


d D M (kpc) 


i q+0-4 


Hi Distance c 


dm (kpc) 
fj, (masyr -1 ) 


1.4 - 7.7 


Proper Motion d 


29 ±1 


Position Angle d 


en 


269 ± 1 



(a) Arzoumanian et al. (1994| ). (b) Taylor fc Cor des ( |1993| ) 
(c) Koribalski et al. ([1995). (d) Bailes et al. ([l990[). 



During the first two nights of this programme (January 
4 and 5, 2001), we discovered an optical bow-shock nebula 
around the radio pulsar B0740— 28. Discovered at Bologna 



(Salter 1970), this pulsar was subsequently shown to be 
rapidly spinning down (McCulloch et al. 1973). Its cor- 



respondingly high spin-down luminosity, E — AirlP/ P 3 , 
(where / = 10 45 /45 g cm 2 is the pulsar moment of iner- 
tia; see Table Q), but likely moderate distance and high 
transverse velocity make it a promising target for power- 
ing an observable Ha nebula. We note that the maximum 
distance derived from H I measurements given in Table | 
is large, but in Section 5 we argue for a distance closer to 
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Table 2. Log of Observations for the PSR B0740-28 field. 



Date 

Instrument 
Pixel Scale ("pix" 1 
Filter 6 (nm) 
Exposure Time (s) 
Transparency 
Mean Seeing (") 
Continuum Frame 



4, 5 January 2001 14 April 2001 



SUSF2/NTT 

0.161 

656/7 

6480 
photometric 

~ 0.9 

None 



TTF a /AAT 

0.373 

656/2 

2700 
photometric 

~ 1.4 
600 s in Rl c 



(a) using the blue Fabry-Perot and R0 (668/21) filter. 

(b) central wavelength and FWHM in nanometers. 

(c) TTF Rl filter (707/26). 



that derive d from the dispersion measure. See Chatterjee 
& Cordes ( 2002 ) for a recent review. 



2. Observations and Reduction 

The nebula around PSR B0740— 28 was discovered in nar- 
rowband frames taken through the 656/7 nm filter of 
SUSI-2 at the 3.5 m New Technology Telescope (NTT), La 
Silla (Fig. Qct). A separate continuum-subtracted narrow- 
band frame (Fig. |^>) was derived from subsequent imaging 
with the Taurus Tunable Filter (TTF; Bland-Hawthorn 
& Jones 1998) at the 3.9 m Anglo-Australian Telescope 
(A AT). Figure 1(c) shows the NTT image in Fig. 0(a) 
smoothed with a 5 x 5-pixel boxcar kernel. The observa- 
tions are summarised in Table |^. Raw frames were treated 
using standard techniques for bias and flat-field correction, 
image registration and co-addition. CCD defects and cos- 
mic rays were removed using cross-pixel interpolation on 
the final frame, in order to preserve as much of the faint 
nebula signal as possible. 

An astrometric solution was derived for the higher- 
resolution NTT image using the iraf task ccmap. A 
Legendre polynomial with 11 coefficients was found to give 
the best fit in x and y to 24 USNO Catalog reference stars 
across the field. The fitting option used in CCMAP sees a 
linear term (with coefficients dependent on translation, 
scaling, rotation and skew), computed separately from a 
distortion term, which is a polynomial fit to the residuals 
of the linear term. No distortion residuals were evident 
in a test transformation of the USNO star coordinates 
to those of the image. The RMS scatter was (Aa, AS) = 
(±0.028",±0.038") and uniform across the field. 

The location of PSR B0740-28 was determined 
by applying the proper motion determined by Bailes 
et al. (1990) to the position measured by Fomalont 
et al. (1992). This places the pulsar at (a, 6) = 
(07:42:49.041, -28:22:44.02) (J2000) at the mean time of 
observation, 1'.'3 inside the leading edge of the nebula 
on its long axis (Fig Q). It is not coincident with the 
prominent star near the head of the nebula. Bailes et al. 
(1990) determine a proper motion position angle (mea- 
sured through east) of — 269 ±1°, i.e. due west. This is 



consistent with the direction of the motion implied by the 
symmetry axis of the nebula and confirms the association. 

3. Flux and Morphology 

The optical emission from pulsar wind nebulae is purely 
Balmer. The Ha flux density in the NTT image was cal- 
ibrated using observations of the planetary nebula flux 
standards 235.3-03.9 and 278.6-06.7 (Dopita & Hua 



1997). Atmospheric extinction was estimated using the 



CTIO measurements of Stone & Baldwin (1983). For the 
planetary nebula standards, the 656/7 filter admits flux 
from the adjacent [N n] A6548 and [N II] A6584 lines, in ad- 
dition to Ha. The [N n] contribution for 235.3-03.9 (37%) 
is more significant than 278.6—06.7 (15%), although the 
656/7 fluxes were corrected for it in both cases. 

The Ha flux over the entire PSR B0740-28 nebula 
is 5 x 10~ 5 photons s _1 cm -2 , around two orders of mag- 
nitude fai nter t han the nebulae powered by J0437 — 471 5 
(Bell et al. |1995[ ) and B1957+20 (Kulkarni & Hester |1988|) . 
The Ha images in Fig. [l] show a distinctive key-hole 
shape which can be divided into three main regions. The 
most western component, or "head" , has an almost cir- 
cular shape, with the pulsar located at the leading edge. 
Approximately 20" from the front of the shock, the nebula 
broadens into a fan-shaped tail and a bridge of emission 
connecting the two sides of the nebula can be seen. At the 
eastern limit of emission, some 45" from the apex, the neb- 
ula appears to show a boundary with opposite concavity 
to that of the apex, closing the tail. Figure 1(6) shows no 
emission within the interior of the head although deeper 
images are needed to confirm this. 

Variations in the surface brightness across the extent 
of the nebula are seen in Fig. 1(a) and Fig. 1(6); see Table 
U for a summary. The knot at Region A, located at the 
juncture of the "fan" and "head" section, and its southern 
counterpart both have surface brightness greater than that 
of the overall nebula. 

4. Nebular Geometry 

When a star moves at a supersonic velocity with respect 
to the interstellar medium (ISM), a bow-shock forms at 
the interface between the stellar wind and the ambient 
medium. For pulsars the wind is relativistic. Assuming an 
isotropic wind, the stand-off distance, r w , separating the 
pulsar from the contact discontinuity of the shock, can be 
derived from the balance between ram pressure and the 
outflowing wind pressure, 



E 



1/2 



471-cpow 2 



(1) 



where po = 2.2 x 10 24 nngcm 3 is the ISM density, rin 



(cm 



is the 



is the hydrogen number density and v p 
space velocity of the pulsar. 

The shock apex was determined by finding the peak 
emission lying along the pulsar proper motion direction 
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Fig. 1. (a) SUSI-2/NTT discovery image (Ha + continuum) of the nebula associated with PSR B0740-28. The arrow 
indicates the direction and distance travelled by the pulsar over 500 yrs. Coordinate epoch is J2000. (b) Continuum- 
subtracted Ha image of the same field taken with TTF on the AAT. Some residuals due to saturated stars and in-focus 
ghost images are present, (c) Image (a) smoothed using a 5 x 5-pixel boxcar kernel. Panels a — c are 1 arcmin on a 
side with north up, east left, with the pulsar position indicated by a circle, (d) Illustrative example of a simple surface 
of rotation based on Eqn. (||) , that assumes the pulsar is moving 20° to the plane of the sky. 



and is located 1'.'3 from the pulsar. However this does 
not correspond to the Balmer emission comes 

from a thin shell upstream from the contact discontinuity. 



Aldcroft et al. ( pL992[ ) find for B1957+20 that the measured 
pulsar-apex separation is 30% larger than r w . Hence we 
adopt r w =1'.'0 = 0.005<ikpc pc, where dk pc is the pulsar 
distance in kpc. The proper motion, when corrected for 
galactic rotation, corresponds to a projected space veloc- 
ity of w p = 1.3 x 10 7 dkp C /cosi cms -1 , where i is the angle 
of this velocity vector to the plane of the sky. Substituting 
v p into Eqn. (|lj) and rearranging gives: 



4cos i L 



n H 



d 4 

"kpc 



45 _3 

cm . 



(2) 



Using the upper and lower limits to the distance 
given m Table g and assuming that i is small then 
0.002 ^ riH ^ 1 cm -3 . If we assume the dispersion mea- 
sure distance c^pc = 2, (which we argue for in Section 5), 
then the density is nn ~ 0.25 cm -3 . 



Table 3. Mean Surface Brightness of Different Regions 



Region 


Surface Brightness 


(see Fig. la) 


(ergs -1 cm -2 arcsec -2 ) 


A 


7xl0 -ia 


B 


4xl0 -19 


C 


5xl0 -19 


entire nebula 


3xl0 -19 



Recent derivations of the geometry of thin momentum- 
balance bow-shocks give a generic form for t he sh ock that 
scales with stand-off distance (Chen et al. 1996; Wilkin 
|1996D - We have compared the shape of the PSRB0740-28 
nebula with the expression of Chen et al. (1996), 



1 



3 

To 



y 



280 r w 



(3) 



Here, our axes lie in the plane of the sky, such that the 
y-axis lies along the proper motion direction, x is perpen- 
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dicular to it and the origin is located at the pulsar. The 
best fit to 14 locations around the "head" gives a stand- 
off distance of l'.'2±0.1 and proper motion position angle 
of 276 ± 3° which are consistent with the measured val- 
ues. However some deviation from an exact fit occur on 
the southern edge of the shock front. Modelling by Wilkin 
(2000) shows that either wind asymmetry or density gra- 
dients in the ISM are able to reproduce such variations. 
Observations of PWNe around other pulsars have demon- 
strated marked anisotropies in the wind flow away from 



the pulsar, (Hester 1998; Gaensler et al. 2002), while large 
variations in the ISM density arc seen in the region of the 



Guitar nebula (Chatterjee & Cordes 2002). 



5. Discussion 

Cordes (1996) shows that the Ha flux density in the nebula 
can be expressed in terms of observables as, 



F n 



10 



-2.47 



Xn H v$8 2 s . 
— — 1 — photons cm s 



(4) 



where X is the neutral fraction, vj — i; p /10 7 and 9 S is 
the observed angular distance. The 2 kpc distance derived 
from the dispersion measure (Taylor & Cordes 1993) is 
consistent with the amount of reddening and closed shape 
of the nebula, in the following ways. 

The observed reddening in the direction of PSR 
B0740— 28 is A v w 1.5 magnitudes at d ~ 2 kpc and dis- 
tances beyond, corresponding to a reduction in Hcv flux 
by a factor of ~ 3 (Neckel & Klare 198C). Combining this 



with a neutral fraction X = 0.01, (reasonable for a warm 
ionised ISM with the nn ~ 0.25 cm~ 3 calculated earlier 
from Eqn. ^ with d^ pc — 2), gives good agreement with 
the observed flux. At larger distances the density derived 
from Eqn. (^|) implies the pulsar is in the hot phase of 
the ISM where the neutral fraction is very small, thereby 
further reducing the predicted flux. 

The direction of the nebula with respect to the plane of 
the sky will also affect the measured flux. Using the shock 
morphology formula given in Eqn. ([^J), we generated a 
three-dimensional model of the shock and then determined 
its shape projected on the plane of the sky for various 
values of i. We find that to best match the observed stand- 
off distance and opening angle i must be small ( ^ 25°) 
and thus its contribution to the low flux is small. 

The overall morphology of the PSR B0740-28 nebula 
is very similar to the Guitar nebula where the "neck" and 



"body" regions (Cordes et al. 1993) are analogous to our 
"head" and "fan" regions respectively. However, there are 
some differences; the head is more circular than the cylin- 
drical neck of the Guitar nebula and the fan is more rec- 
tilinear than the bulbous body of the Guitar nebula. We 
note, however, that both nebulae show stronger Hcv emis- 
sion in the transition between the narrow and broad re- 
gions. These could indicate regions of increased density on 
sub-parsec scales, given that the full extent of the nebula 



is ~ 0.6 pc if 2 kpc away. Alternatively, the stronger emis- 
sion may indicate the transition between different physical 
regimes of the nebula. 

The Guitar nebula emission shows two distinct bub- 



ble regions (Cordes et al. 1993). The fainter, inner bubble 
is analogous to the "bridge" of emission that we see in 
PSR B0740—28. It is possible to generate these apparently 
closed structures via projection of a truly open structure. 
This can occur if the nebula is slightly inclined, as is shown 
in Fig 1(d) for a cone assumed to be inclined at 20°. 

The apparently closed nature of the PSR B0740-28 
nebula suggests that in the rear region of the nebula, the 
ram pressure due to the pulsar motion is no longer im- 
portant and we could instead be seeing the point where 
the pressure in the wind is balanced by the ISM pres- 
sure — it is possible that the break between the "head" 
and "tail" components of the PWN may correspond to 
the transition between these two sources of confinement. 
However, this explanation is not consistent with the en- 
tire PWN being embedded in gas of uniform density, since 
the pulsar's motion must be supersonic to form a bow 
shock, but in which case there can be no pressure bal- 
ance at the rear of the PWN. In order to explain the 
observed morphology, we propose that the pulsar has re- 
cently passed through a strong density gradient in the 
ISM, such that the pulsar's motion is supersonic in the 
comparatively dense medium (tih ~ 0.25 cm~ 3 ) through 
which it is currently propagating, but was subsonic in 
low density material (nn 0.003 cm~ 3 ) surrounding the 
closed end of the PWN. Similar arguments can be invoked 
to explain the morphology of the Guitar PWN powered by 
PSR B2224+65. 

Considering the latter region, the requirement that the 
pulsar be moving slower that the local sound speed im- 
plies an upper limit to the pulsar's distance of 3-4 kpc, 
as no ISM component has a sound speed faster than 
- 500 km s _1 . If PSR B0740-28 is at 2 kpc, then the 
pulsar wind pressure, P w = E / (Airr^c) approximately bal- 
ances the ISM pressure, Pism ~ 5 x 10~ 14 Pa, at the 
the location of the closed region (where is the distance 
to the back of the nebula), providing further support for 
the lower distance argued for above. In the interpretation 
we propose, the closed region of the PWN corresponds 
to the static or "ghost" nebulae proposed by Blandford 
et al. ( 1973 ). The faintness of this region is consistent 
with the expectation that such sources should generally 
never be observable (Gaensler et al. 200C); in the case 
of PSR B0740-28, it was only through the presence of 
brighter emission near the pulsar that this closed region 
was identified. 



Bucciantini & Bandiera (2001) have recently shown 
that classical thin-shock models are not appropriate for 
pulsar bow-shock nebulae. However, here we find reason- 
able agreement between the classical models and the shape 
of the head of the PSR B0740-28 nebula. It is likely 
though that the microphysical processes they discuss are 
responsible for at least some of the morphological charac- 
teristics of these PWNe. Indeed Bucciantini & Bandiera 
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( |200lD point out that perhaps PSR B0740-28 is an in- 
termediate case between the nebulae associated with the 
MSPs and the Guitar nebula. Interestingly, they also pre- 
dict that, based on the ISM conditions, this nebula should 
be readily detectable. However, we find that it is con- 
siderably less-so than the others. Clearly, increasing the 
presently small sample of optical bow-shock nebulae re- 
mains a key step towards being able to generalise about 
the ISM conditions within the vicinity of pulsars. 
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